Stripes in thin ferromagnetic films with out-of-plane anisotropy 
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We examine the T — phase diagram of a thin ferromagnetic film with a strong out-of-plane 
anisotropy (e.g. Co/Pt multilayers) in the vicinity of the reorientation phase transition. The phase 
diagram in the anisotropy-applied field plane is universal in the limit in which the film thickness is 
the shortest length scale. It contains uniform fully magnetized and canted phases, as well as period- 
ically nonuniform states: weakly modulated spin-density waves and strongly modulated stripes. We 
determine the boundaries of metastability of these phases and point out the existence of a critical 
point at which the difference between the SDW and striped phases vanishes. Out-of-plane magneti- 
zation curves exhibit hysteresis loops caused by the coexistence of one or more phases. Additionally, 
we study the effect of a system edge on the orientation of nearby stripes. We compare our results 
with recent experiments. 
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I. INTRODUCTION 

Recent experimental studies have revealed rich physics 
of thin ferromagnetic films with an easy axis of magne- 
tization normal to the film plane»i>2»2»£iiii£»ii& Such films 
possess several magnetic phases in which the magneti- 
zation can be uniformly normal to the plane, canted, or 
periodically modulated in one direction (striped). Mag- 
netization curves exhibit intricate magnetic hysteresis in- 
dicating coexistence of phases. Magnetic probes with 
sub- micron resolution provide detailed information about 
nucleation and growth of new domains during the mag- 
netization reversal^ These developments fuel the need 
for a theoretical understanding of ferromagnetism in thin 
films. 

Theoretical studies of thin films with magnetic mo- 
ments pointing out of the plane date back to the 1980s. 9 
The simplest model describes Ising spins with local ferro- 
magnetic exchange and long-range dipolar interactions. 
A competition between these forces makes the uniform 
ferromagnetic states unstable towards a spontaneous for- 
mation of magnetic stripes with an alternating sign of 
the magnetization. The stripe period is a mesoscopic 
length determined by the relative strengths of exchange 
and dipolar interactions. 

A major drawback of the dipolar Ising model is the 
neglect of the in-plane components of magnetization, 
which become important when the magnetization ro- 
tates away from the plane normal, a process known as 
the reorientation phase transition (RPT). A minimal 
model must therefore use a three-dimensional vector of 
magnetization whose magnitude M is considered fixed 
well below the Curie temperature and whose orienta- 
tion is given by the colatitude 9 and azimuth <f>: M = 
M(srrj,0cos</>, sin 9 sin 0, cos#) with the z axis normal to 
the film plane xy. In a state with uniform magnetiza- 
tion, the energy density contains the easy-axis anisotropy 
—Kcos 2 9 and the demagnetizing term (/io/2)M 2 cos 2 9 
due to the magnetic field. In the absence of an applied 
field the RPT occurs when the anisotropy drops below 
the critical strength Kq = [IqM 2 /2. 



The assumption of uniform magnetization breaks down 
in the vicinity of the RPT: as in the dipolar Ising model, 
the competition between local and long-range forces leads 
to the formation of stripes with a mesoscopic period that 
depends on the effective anisotropy K — Kq, exchange 
strength A, and the film thickness t^ Depending on the 
anisotropy and on the strength of an applied magnetic 
field, the stripes can appear as a weak spin-density wave 
(SDW) in the background of uniform magnetization or as 
fairly wide domains of uniform magnetization separated 
by narrow domain walls. 

The presence of the in-plane components of magneti- 
zation leads to an important distinction of domain walls 
from their counterparts in the Ising model: the vector of 
magnetization rotates between the (mostly) upward and 
downward directions. Thus domain walls are endowed 
with in-plane magnetization, a fact with important topo- 
logical consequences. The domain walls are of the Bloch 
type: the in-plane magnetization on them points along 
the wall. In contrast to Neel walls (in-plane magneti- 
zation normal to the wall), Bloch walls do not generate 
stray magnetic field and thus have a lower magnetic en- 
ergy. We will therefore specialize to magnetization con- 
figurations in which the vector of magnetization depends 
on a single coordinate x and lies in the yz plane: 



M = M(0, sin0(x), cos0(x)). 



(1) 



Using a variational approach, Berger and Erickson 11 
obtained a phase diagram of such a one-dimensional 
model as a function of the anisotropy K — Kq and an 
applied out-of-plane field H±. It exhibits several phases 
with both first and second-order transitions between 
them. Berger and Erickson focussed on thermodynamic 
transitions and did not provide boundaries of metastabil- 
ity. Such boundaries are important for the understanding 
of magnetic hysteresis in thin films. Often thermal acti- 
vation is insufficient to initiate the decay of a metastable 
phase and magnetization reversal begins only when that 
phase becomes locally unstable. 

In this work we describe several new results. Our main 
achievement is the determination of out-of-plane magne- 
tization curves M±(H±) that can be directly compared 
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with experimental data. We show that, for sufficiently 
thin films, the shape of the magnetization curve depends 
on a single parameter that is a function of the dimension- 
less effective anisotropy k = (K — Kq)/Kq, film thickness 
t, and the magnetic exchange length A = \/2A/ '[1qM 2 . 
This universality is the consequence of a scaling prop- 
erty of the free energy in the thin- film limit. A proper 
rescaling of the anisotropy k and magnetic field H± yields 
a universal phase diagram. We point out the existence 
of critical points of the liquid-gas type that terminate 
lines of thermodynamic first-order phase transitions be- 
tween striped and SDW phases. Finally, we discuss the 
behavior of stripes near the edge of the film. We find 
a tendency for stripes to meet an edge perpendicularly, 
independent of the particular shape of stripe. 

The formalism used in this work is minimization of the 
free energy. The average out-of-plane magnetization in a 
sample with stripes is chiefly dependent on the average 
period of those stripes rather than their orientational or 
translational order. We thus neglect the influence of ther- 
mal fluctuations, which tend to disrupt that orientational 
and translational order in the stripe phase i 12 ' 13 

The paper is organized as follows. In Section [TT] we 
derive the functional of magnetic free energy special- 
ized to one-dimensional variations, describe its scaling 
properties in a thin-film limit, and introduce appropri- 
ately rescaled variables. We illustrate its use by find- 
ing the lines of instability of the uniform phases in the 
anisotropy-applied field phase diagram. In Section [TTT1 we 
determine the boundaries of metastability of nonuniform 
phases and magnetization curves using a numerical min- 
imization. We also find the location of the stripe-SDW 
critical points in the (K, H±) plane. Section [TVl deals 
with the behavior of stripes near the film edge. Some use- 
ful intermediate results are described in the Appendixes. 



where the dipolar kernel V^r) = 1/r — l/\/r 2 + f 2 re- 
flects the interaction of magnetic charges with densities 
±Mm z (r) induced on the top and bottom surfaces of the 
film. The expression for the dipolar energy ([3]) is valid 
provided that there are no magnetic charges in the bulk 
of the film, i.e. d x m x + d y m y = 0. This condition is 
compatible with Eq. |T]), which describes a system with 
Bloch domain walls. Domain walls of the Neel type gen- 
erate additional dipolar terms. 

Specializing to one-dimensional configurations without 
bulk magnetic charges ([I]), and with magnetic field ap- 
plied perpendicular to the plane of the film, we obtain 
the energy 
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A(d9/dx) - Kcos 2 9- [i MH cos 6) 



+ -fi M 2 j dx j dx' cose(x)V{x-x') cose(x'),(A) 

where L y is the width of the system in the y 
direction. The one-dimensional kernel V{x) = 
(l/27rf)ln(l + t 2 /x 2 ) obtained by integrating (O over y 
has the Fourier transform 
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The Taylor expansion (0 is justified when the film thick- 
ness t is the shortest length scale in the problem. The 
zeroth-order term Vo(k) — 1 can be interpreted as a con- 
tact part of the dipolar interaction that simply shifts 
the anisotropy: K >—> K — Ko- The first-order term 
Vi(fc) = —\k\t/2 represents the effect of the stray dipo- 
lar field. Its inverse Fourier transform V\{x) diverges at 
short length scales and requires a short-range cutoff. It 
has the following properties: 



II. THE MODEL 
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oo, / dxVi(x) = 0. (6) 



A. The free energy 

The free-energy functional for a thin ferromagnetic film 
of thickness t can be separated into a local and long-range 
parts. The local part includes the exchange, uniaxial 
anisotropy, and Zeeman energies: 



-Blocal — t 



d 2 r (A\ Vrh| 2 - Krn 2 z - fj, M H • m) , (2) 



where m = (sin 8 cos 4>, sin 9 sin <fi, cos 6) is the three- 
dimensional unit vector pointing along the magnetization 
and V = (d x , d y ) is the two-dimensional gradient in the 
plane of the film. The long-range part is due to dipolar 
interactions: 
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J d 2 r J dV m z (r)V(r-r')m z (r'), 



(3) 



B. Scaling property of the free energy 

It is convenient to use natural scales for various phys- 
ical units: Ko — noM 2 /2 for the effective anisotropy k 
and volume energy density, and M for magnetic field: 



k= (K-K Q )/K , 
In these units, 



H/M. 



(7) 
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dx' cos #(2;) Vi(x — x') cos#(x'), 

(9) 
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where A = yj A/Kq is the exchange length. 

The free energy, given by the sum of these terms, is 
invariant under a scaling transformation 

x i ► bx, t h-> t, An 6 1 / 2 A, k i ^ h i— > 6 _1 h. 

(10) 

This symmetry indicates that the state of the film de- 
pends on the effective anisotropy and magnetic field 
through scale-invariant variables k/kq and h//co, where 
Ko = f 2 /(4A) 2 is an effective anisotropy scale whose sig- 
nificance will be clarified shortly, ft also yields a charac- 
teristic length scale 8tt\ 2 /t, which determines the period 
of the stripes — see Eq. (|12p below. 

The free energy is scale- invariant only to the lowest or- 
der in the film thickness t. Inclusion of higher-order terms 
in the dipolar kernel §5§ violates this property. The scal- 
ing works as long as the thickness t is small compared to 
all other length scales, in particular the exchange length 
A«5 nm (for Co). 

Note that this scaling law applies more generally to 
two-dimensional configurations with an applied field in 
any direction, as the energies associated with bulk mag- 
netic charges and in plane fields are also invariant under 
flTUJl when t is small. 

C. Uniform phases 

Let us first discuss the uniform phases and their insta- 
bilities in the case that the applied field is normal to the 
plane. The case k > is trivial: only the upward and 
downward polarized states are stable. Both are locally 
stable in the region \h\ < k. The coexistence leads to a 
hysteresis in magnetization curves m z (h). 

The situation is more interesting for k < 0, where 
magnetization prefers the in-plane direction. In a uni- 
form state the stray field vanishes and the energy is pro- 
portional to — (k/2)cos 2 6* — hcosO. In a strong field, 
\h\ > \k\, the film is fully polarized, cos 6 — sgn(/i). Be- 
low the critical strength, \h\ < \k\, the magnetization is 
canted, cos# = h/\n\. The uniform-to-canted transitions 
at h = ±K are continuous. 

Consider the free energy of small fluctuations around 
a canted state, cos8(x) — —h/n + S(x): 

1 J Zir \sm I ) 

(11) 

where 5(k) is the Fourier transform of 8 (x) and cos#o = 
—h/n is the equilibrium value for the canted state. The 
softest mode has the wavenumber 

t sin 2 6>p 

k ° = 4A~2 ' (12) 

A spin-density wave develops in the canted background 
on the line 

{h/nf = 1 + «/kq, (13) 



where Ko = £ 2 /(4A) 2 . An expansion to the order 0(S 4 ) 
reveals a positive-definite quartic term. Thus the canted- 
SDW transition is also continuous. Note that the scaling 
arguments are confirmed: the wavelength is indeed set 
by the scale 8ttX 2 /t and the critical line (fT3]) contains the 
rescaled variables h/ k$ and k/kq. 
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FIG. 1: (Color online) Phase diagram showing the regions of 
stability of uniform phases. Up (U) and Down (D) phases 
coexist in the region on the right. Solid and dashed lines 
denote continuous and discontinuous transitions, respectively. 

Since no uniform solution is stable inside the semicu- 
bic parabola (JT3J) , this region must be occupied by an 
inhomogeneous state, which one might easily guess to 
be a stripe phase (see Fig. [TJ. However, the situation 
is somewhat more complicated. A further analysis shows 
that the stripe phase remains (at least locally) stable out- 
side the semicubic curve. In addition, we find a region 
of coexistence between a strongly inhomogeneous striped 
state and a weakly inhomogeneous SDW state. Like a 
liquid and a gas, the two phases differ from each other 
)) only quantitatively. Indeed, there is a critical point at 
which the differences vanish continuously, as we will dis- 
cuss below. 



III. STRIPE PHASE 

In this section, we describe the striped state in various 
regions of the h—K plane. We determine the approximate 
boundaries of stability of the striped phase, and show 
how the metastability of the stripe phase leads to the 
hysteresis curves observed in experiments. 
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Numerical determination of metastability limits 
for the stripe phase 
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FIG. 2: (Color online) Numerically determined boundaries of 
the nonuniform phases are shown as symbols. The Striped 
phase can coexist with the Canted (C), Up (U), Down (D), 
and Spin density wave (SDW) states in various regions of the 
phase diagram. Solid and dashed lines denote continuous and 
discontinuous transitions, respectively. 

In order to find the boundaries of stability of the 
striped phase, we performed a numerical simulation of 
the system described by Eqs. (|8|) and ([9]). The simu- 
lation was conducted on a chain of 2048 magnetization 
vectors of unit length described by the angle 9(x). Peri- 
odic boundary conditions were applied. We used specific 
values for the exchange length and thickness of 5 and 3 
lattice spacings, respectively. However, based on the scal- 
ing arguments described in Section III B[ we expect the 
results to be universal for thin films. Eq. ([5]) was mini- 



mized using a relaxation method ^ where h was replaced 
by h e g(x) = h— J dx' Vi(x — x') cos9(x') to account for 
the dipolar stray field. This h c g was recalculated after 
a number of iterations of the relaxation method. The 
process continued until the Lagrange equation, 

A 2 (d 2 6»/dx 2 ) = rccos0sin0 + /i e ffsin0, (14) 

was satisfied at each point. Random noise was then 
added to the system and h was incremented. In this way, 
the algorithm moved the state of the system along the 
local minimum of energy. Sweeps of the magnetization 
were conducted from positive to negative values of h and 
the average magnetization was recorded. The hysteresis 
loops shown in Fig. [3] were produced from these single 
sweeps by rotating the data points to produce the up- 
ward sweep and superimposing it on the downward one. 
Discontinuous jumps in the magnetization mark the sys- 
tem's entry into and exit from the striped phase. In what 
follows, we interpret the resulting phase diagram (Fig. [5]) 
for various values of the effective anisotropy. 

An important feature of phase diagram revealed by the 
numerical simulation is the existence of liquid-gas like 
critical points at k c w — 0.82ko, h c w ±0.3ko (Fig. 
To the right of this point, a distinction can be made 
between a striped phase, with large modulation around 
a small average magnetization, and a spin-density wave 
phase, which has a small modulation around a larger av- 
erage magnetization. Boundaries of metastability extend 
from the critical points, surrounding regions of coexis- 
tence between the striped and SDW phases. To the left 
of k c the two phases merge into one and magnetization 
curves proceed in a reversible manner. It is important to 
note that for n c < n < 0, the second-order phase tran- 
sition line out of the canted phase is very close to the 
first-order phase transition line into the striped phase. 
As such, this second-order transition may be difficult to 
detect as it is hidden by the nearby first-order jump. 

In the remainder of this section, we describe the results 
of a sweep of the applied magnetic field beginning in high 
positive field normal to the sample and moving to a high 
negative field. The results of such a sweep depend on the 
anisotropy constant k/kq describing the system. 

For strong enough in-plane anisotropy, n < —kq, the 
striped phase is never seen. A sweep of magnetic field 
from high positive to high negative fields would find a 
completely reversible magnetization curve, with the uni- 
form upward phase beginning to cant at h = — k and 
following the field smoothly and linearly. The magneti- 
zation will be —h/n until the field is decreased to the 
h = k line, from which point onward the sample is down- 
wardly polarized. 

For anisotropy values only slightly greater than the 
RPT value n = — Kq, a modulation of the magnetization 
will appear in low fields. The magnetization curve will 
still be reversible, but no longer exactly li near. Onc e the 
magnetic field crosses the curve h — — fcyl + k/kq from 
above, a small modulation out of the canted phase de- 
velops continuously. That modulation fades again at the 
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lower curve h = «yl + k/kq, and the system proceeds 
as before from the canted to the uniform down phase. 
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FIG. 3: (Color online) Numerical simulation of the out-of- 
plane hysteresis curves for several values of k/kq. The lo- 
cations of these sweeps in the phase diagram are shown by 
vertical lines in the bottom panel. 

For anisotropy values greater than k c , the system 
shows history dependence in the magnetization curve. 
At a slightly smaller value of h than h — —Ky/l + k/kq, 
while the amplitude of the spin density wave is still on the 
order of 10 -3 , the system undergoes a first-order transi- 
tion from a state with small modulation and large average 
magnetization to one with small average magnetization 
and a large amplitude modulation. This stripe phase 
remains as the field decreases until finally the system 
undergoes a first-order transition in which the periodic 
modulation disappears. The boundary of stability of the 
striped phase crosses h = k, the line beyond which the 
downwardly magnetized state is stable, so that for more 



negative values of k, the striped phase decays to a canted 
phase, while for more positive K values, the system enters 
the downwardly polarized state. For example, in Fig. [3j 
the k — — 0.81ko and k — — 0.65ko magnetization curves 
show the striped phase decaying into the canted phase, 
while the other loops show a decay directly into the fully 
polarized phase. 

For a small k > 0, as we conduct a downward sweep 
of the field beginning at large positive values, the sys- 
tem undergoes a first-order phase transition out of the 
upwardly polarized state at h = —k. There is no canted 
or spin-wave phase after the transition, but the simula- 
tions indicate that there is an intervening stable state 
with non-constant magnetization before the system de- 
cays into the downwardly polarized state. The limita- 
tions of the simulation become evident here. The peri- 
odic boundary conditions imposed in the simulation do 
not allow us to distinguish between a periodic state with 
a long period and a state with an isolated region of un- 
favored magnetization (soliton). Only a small number of 
solitons appeared within the period forced by the simu- 
lation. The considerations of Sections IIIIB H and If 11 B 2l 
indicate that for large enough k the isolated soliton state 
is stable while the striped state is not. Whether this 
holds at small k > has not been determined. 

The non-uniform state resulting from the first-order 
phase transition out of the upwardly polarized state per- 
sists until it vanishes in a first-order transition to the 
downward state (at h « — 0.5ko m the simulation). If 
one were to reverse course and increase the field again be- 
fore reaching this second transition, the upward regions 
that had been unfavored would grow slowly rather than 
the system returning immediately to the totally upwardly 
polarized phase. This behavior accounts for the "fading 
contrast" seen in experiments with Co/Pt multilayers^ 

If ft > is large enough, we expect that the transition 
will proceed directly from an upwardly polarized state 
to a downwardly polarized state at h = —k. For very 
large k, neither a striped nor an isolated-soliton phase 
can form from the decay of the upwardly polarized state, 
as neither will be stable in the high fields required for the 
uniform state to become unstable. (See Sections IIII B II 
and|IIIB2J). 

There is a caveat on this expectation, however. If topo- 
logically non-trivial domains of reversed magnetization 
form during the transition out of the upwardly polar- 
ized phase, it may be that they are stable even in these 
large fields. A topologically non-trivial domain in our 
case is one for which the magnetization angle 9 at the 
left end of the domain differs by a non-zero multiple of 
2n from the magnetization angle at the right end. It 
is not unreasonable to expect the formation of such do- 
mains: if two downward regions form during the transi- 
tion by rotating in opposite directions away from 8 = 
then the upwardly polarized region between them will 
be topologically non-trivial. Since these downward re- 
gions may be initially well separated, their rotation di- 
rections are essentially independent, and so a non-trivial 
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upward domain will form between them roughly half the 
time. These non-trivial domains are the one-dimensional 
analog of skyrmion-type domains in two dimensions. 15 
Skyrmions may be the cause of the asymmetric domain 
nucleation observed in Co/Pt multilayers 

The magnetization curve would look nearly the same 
in the case that topologically non-trivial solitons form 
as it would if they do not. The topologically non-trivial 
solitons would have small width when the field is large 
and so change the magnetization by a negligible amount. 
However, when the field is reversed after the transition, 
these soliton domains would grow. This would cause the 
system to proceed through the striped phase rather than 
staying downwardly polarized until h — k. Hysteresis, 
then, would not be observed unless the applied field is 
strong enough to force the width of topological solitons 
down to the lattice scale, allowing them to decay. 



B. K>Ko, the wide-stripe approximation 

When the out-of-plane anisotropy is strong, the mag- 
netization pattern is expected to consist of long up- or 
downwardly pointing regions, separated by short regions 
in which the magnetization changes rapidly from one do- 
main type to the other. We refer to such a configu- 
ration as a "wide-stripe" phase. (Note, however, that 
this "phase" is continuously connected to the modulated 
phase at negative k, and so not a truly distinct phase.) 

It is well known that an energy advantage can be 
gained over a purely uniform phase through such modu- 
lation. In fact, when the out-of-plane anisotropy is very 
strong, the zero field model differs little from the dipo- 
lar Ising model discussed by Garel and Doniach. 9 The 
two significant differences lie in the energy costs of the 
kinks and in topological considerations. The topological 
consequences of the model will be discussed in Section 

hub 21 

The wide-stripe configuration can be thought of as a 
periodic array of Bloch domain walls. These walls can 
be treated as elementary objects with some internal en- 
ergy cost and long range interactions with other walls 
through the stray dipolar field. The interaction energy 
of two walls is logarithmic in their separation. If we de- 
fine the orientation of a wall as a vector pointing from the 
upwardly magnetized side of the wall to the downwardly 
magnetized side, we find that the magnetic interaction 
is attractive for walls with the same orientation and re- 
pulsive for walls of opposite orientation. For a periodic 
structure, walls always have orientation opposite that of 
their nearest neighbors. It is this nearest neighbor repul- 
sion caused by the stray magnetic field that allows for 
the stability of a striped phase despite the cost of the 
walls in exchange energy and the exchange force acting 
between nearby walls. 

These considerations can be used to derive an analytic 
expression for the magnetization curve in the wide-stripe 



phase: 
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where we have used the rescaled variables 

n ~ h 
k = — , h = — . 
n k 



(15) 



(16) 



Note that this magnetization curve is consistent with the 
scaling property described in Section III Bl Details of the 
derivation, including an analytic expression for the stripe 
period, can be found in Appendix [A"l 



1. Zero-soliton-density transition 

As the external magnetic field is varied, the wide- 
stripe phase becomes dominated by the regions of mag- 
netization favored by the magnetic field. For a strong 
enough magnetic field the striped phase has only nar- 
row, widely separated regions of the unfavored magne- 
tization (solitons). As Eq. (|A15|) shows, the period of 
the wide-stripe phase tends to infinity along the curve 
h = — (4/7r 2 )\/^exp(— That is, the density of 
the unfavored solitons reaches zero. For fields above this 
curve, there is no stable structure with evenly spaced do- 
mains of unfavored magnetization. Note that the field at 
which the solitons are expelled from the system decreases 
exponentially with k. For large k the region in which 
stripes exist is extremely narrow. As the field varies from 
down to up across this region, a system that is nearly 
all down with a few isolated upward solitons will move 
quickly through the striped phase to a state that is nearly 
all up with a few isolated downward solitons. 



2. The collapse of isolated solitons 

In a perfectly isotropic sample, the field drives all the 
regions of unfavored magnetization to the edge of the sys- 
tem during the zero-soliton-density transition described 
above. A real sample, however, may have pinning cen- 
ters where regions of opposing magnetization would be 
localized. These isolated solitons will persist into a much 
higher field, and their collapse is dependent on their topo- 
logical character. 

The edges of such regions interact with each other 
through two distinct forces. There is a long-range repul- 
sion between them due to the dipolar stray field and a 
short-range exchange interaction that may be attractive 
or repulsive depending on the topology of the region. In 
addition, there is a force from the applied magnetic field 
that acts to increase or decrease the width of the soliton 
based on its polarization. If the exchange force is at- 
tractive, the boundary of stability of the isolated soliton 
phase is found at the field that will squeeze the soliton 
boundaries enough that the exchange force takes over and 
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the solitons collapse. If the exchange force is repulsive, 
the soliton must be forced to a width smaller than the 
lattice spacing in order to collapse, so that a phase dif- 
ference of 2-7T between neighboring sites may be ignored 
as physically meaningless. 

The exchange force between two domain walls a dis- 
tance w apart is 



F ex (w) ~ ±8«;exp 



as 



(17) 



where w — t/kqw/X (see Appendix [Bj. The force is 
attractive if the kinks form a non-topological soliton and 
repulsive if the soliton is topological. 

To determine the field necessary to collapse a non- 
topological soliton, we take the Zeeman energy of the 
soliton and the energy of the stray field to be 



IV. STRIPES NEAR AN EDGE 

In this section, we consider the orientation of stripes 
near a system edge. Although we begin with a variational 
solution with a sinusoidal modulation, the final result is 
applicable to a general profile for the magnetization. We 
find that stripes oriented with their domain walls perpen- 
dicular to the edge are energetically favored over all other 
orientations except possibly that with stripes exactly par- 
allel to the edge. For simplicity, we consider only the case 
of no applied field. The results are entirely similar when 
a field is applied. In particular, there is no change to 
Eq. (f2l)]) . below. We compare these results with recent 
experiments on thermally evaporated Ni films. 

We use the trial solution 



a sin (q • x — (3) . 



(23) 



E(w) = const + 2hw lnw; 

7T 



(18) 



for a single downward soliton of width w [sec Eq. (IA20[) ]. 
We use here the \fkw 3> 1 approximation, as we expect 
the collapse width of a non-topological soliton to be large 
compared to the size of the walls that border the soliton. 
As the field h increases, the soliton is squeezed: 



2h 



8«;exp 



(19) 



Note that as w shrinks, the restoring force on the RHS of 
Eq. (| 19(1 increases at first and then reaches a maximum 
at 



In k 
2 a/I 



as 



(20) 



A further increase in the field leads to a collapse of 
the soliton, as the restoring force can no longer bal- 
ance the force due to the field. This justifies the ap- 
proximation of large width, since for the collapse width 
y/kw ~ (1/2) 1iik> 1. Hence, in a field 



8y/k 
7r In k 



as 



oo, 



(21) 



the (non-topological) soliton will collapse. 

If, however, the soliton consists of a full 2ir rotation 
of the magnetization, the restoring force never reaches a 
maximum. While the dipolar contribution to the force ef- 
fectively disappears as the soliton width decreases to the 
order of the domain wall width, the repulsive exchange 
force increases without bound. The soliton can then only 
collapse when its width 



w= w = 2\ arcginh LM (22) 

y/KO y/K+\h\ VV \h\J 

derived in Appendix IB! reaches the lattice scale. 



Here a is the amplitude of oscillations; the stripe 
wavevector q = (qcosa, gsina) has a fixed length q; 
the angle in the plane a ranges from for stripes normal 
to the edge to ir/2 for stripes parallel to the edge; the 
phase P will be important only for stripes parallel to the 
edge. In order to represent an edge at y = 0, we include 
a a step function <d(y) in all y integrals. 

The energy associated with the presence of an edge 
comes from two sources. First, stripes at the edge gen- 
erate a stray magnetic field. Second, depending on the 
sign of the effective anisotropy k, the system can lower 
the energy by having an extra node or antinode near the 
edge. (The latter works only for stripes parallel to the 
edge.) 

The long-range nature of the dipolar forces requires 
a certain amount of care in evaluating the edge energy. 
The value of the dipolar energy is sensitive to three length 
scales: the film thickness t, the width L y , and the stripe 
period 2n/q. On a computational level, keeping t and 
L y finite is required to avoid ultraviolet and infrared di- 
vergences. Fortunately, the difference in the energies be- 
tween stripes with different orientations (parameterized 
by the angle a) is insensitive to these length scales, as 
long as a ^ ir/2. This simplifies the computation greatly. 
Details of the calculation can be found in Appendix [C] 

For a ^ 7r/ 2, the energy difference 



AE = E(a,q)-E(0,q) 



(24) 



has no ultraviolet or infrared divergences, and is in fact 
independent of q (Fig. [4j) : 

AE(a) 



H M 2 t 2 L x 8tt 



[sin a In (1 + sin a) + (1 — sin a) In cos a]. 

(25) 

If a = tt/2, the stripes are oriented exactly parallel to 
the edge. In that case there is a cutoff dependent term 
in the stray field energy difference proportional to 



8tt 



In (qt) cos 2/3. 



(26) 
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As long as the wavelength of the mode under consider- 
ation is larger than the thickness, the logarithm is neg- 
ative. This term is therefore minimized when /3 = n/2 
and so acts to attract the crests of the spin wave to the 
system edge. 

There is also an additional term associated with the 
local part of the free energy that appears only when a = 
ir/2. It is proportional to 

a 2 /, 9 9 \ sin (2/3) . . 

- T (AV + *)-^. (27) 

If the coefficient in brackets is positive (negative), this 
term acts to minimize the number of nodes (crests) of 
the spin wave that are present in the system, since each 
node (crest) has a cost in the exchange and anisotropy 
energy. The competition between the terms (|2T)1) and (j2"T|) 
will determine the phase of a spin-density wave parallel 
to the edge. 




(27) 



1 



0.8 - 




stripe angle 

FIG. 4: Energy as a function of the angle the stripe wavevec- 
tor makes with the system edge. Energy here is measured in 
units of (a 2 /8n)fj,oM 2 t 2 L x . 

The addition of the phase-dependent terms above may 
mean that stripes parallel to the edge end up being the 
global minimum of energy. However, this minimum takes 
the form of a downward spike in the energy that occurs 
at an angle that would otherwise be a maximum of the 
energy (see Fig.|4|). Even if the state with stripes parallel 
to the edge is the global minimum of energy, a low tem- 
perature system that begins with stripes at a random ori- 
entation is more likely to evolve toward the (metastable) 
minimum at a = 0. 

Since the angular dependence of the energy for a ^= n/2 
is independent of the wavenumber q, this result holds 
for any stripe profile that may be made up of higher 
harmonics. In that case, a 2 /2 should be interpreted as 
the average value (cos 2 9) of the striped state. The result 
is also independent of any of the material parameters 
because it is the energy of the stray field that causes the 
effect. 

This result is consistent with observations made in 
thermally evaporated thin nickel films - and numerical 



FIG. 5: (Color online) A stationary configuration obtained 
from a random initial state in a disk of thickness 14 nm and 
diameter 400 nm. Magnetization length M = 1.4 x 10 6 A/m, 
exchange constant A — 3.3 x 10~ n J/m, exchange length 
A = 5.2 nm, easy-axis anisotropy K — 1.5 x 10 6 J/m 3 yield 
K = 0.22 and = 0.45. Magnetization points up in the red 
(light gray) regions and down in the blue (dark gray) regions. 
Numerical simulation using OOMMF.— 



simulations (Fig. [5]) that indicate that walls meet the 
edges of the film at a 90° angle. 

V. DISCUSSION 

We have examined the properties of thin ferromag- 
netic films with a strong easy-axis anisotropy K w 
Kq = HqM 2 /2 favoring the out-of-plane component of 
the magnetization, with Co/Pt films as a prototype. 
We have found a scaling property that applies in the 
limit where the film thickness is the shortest length 
scale in the problem. In such a case the phase diagram 
in the field-anisotropy plane is universal if the applied 
field H and effective anisotropy K — Kq are properly 
rescaled. The proper variables are h = H/(ko,M) and 
k = (K — Kq)/(k Kq). The parameter k q = (i/4A) 2 
is determined by the film thickness t and the exchange 
length A = y/2A/fi Q M 2 . 

The universal phase diagram in the case of an ap- 
plied field normal to the plane was determined through 
a combination of analytical and numerical methods fo- 
cussed on uniform states and on states with a one- 
dimensional modulation of the magnetization, such as 
magnetic stripes observed near the reorientation phase 
transition (K ~ Kq). In addition to fully magnetized 
and canted uniform states, at least two non-uniform mag- 
netic states were found: a spin-density wave and a striped 
phase. These two states are found to coexist in parts of 
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the phase diagram (like a gas and a liquid). The co- 
existence of various phases (e.g., stripes and SDW or 
stripes and canted) explains the rich hysteretic behav- 
ior of magnetization observed in these films. We have 
determined the boundaries of stability and metastability 
of these phases and obtained out-of-plane magnetization 
curves M±(H±), which have universal shapes determined 
by the rescaled anisotropy k. 

In addition to developing the zero-temperature phase 
diagram near the reorientation phase transition, we have 
expanded on previous work in the case of large out-of- 
plane anisotropy by finding analytic expressions for the 
stripe period and magnetization as functions of the ap- 
plied field. We have shown that the range of field values 
for which a structure of evenly spaced stripes is stable is 
exponentially small for large anisotropy. The period of 
such a structure tends to infinity at a small value of the 
applied field. Stripes arc unlikely to be found, then, at 
large anisotropies. 

We have investigated the behavior of stripes near the 
film edge and found that the energy of the dipolar stray 
field is minimized when the stripes are perpendicular to 
the edge, as recently observed in Ni films. At the same 
time we find that a state with stripes parallel to the edge 
provides an opportunity to lower other terms in the free 
energy, e.g. the magnetic anisotropy, by registering the 
nodes or antinodes of magnetization at the system edge. 

The behavior of the striped phase in an out-of-plane 
magnetic field depends on the topology of the stripes, 
i.e. on the direction of rotation of the magnetization. If 
in the absence of an applied field the colatitude angle 9 
oscillates between (nearly) and it, the topology is triv- 
ial. As the strength of the applied field increases, the 
stripes of the "wrong" orientation first shrink and then 
disappear leaving behind a uniform state. However, if 
the domain walls delineating a "wrong" domain wind in 
the same direction (say from to n and from n to 2tt), 
this domain has a nontrivial topology. As its wall are 
pushed together, they feel strong repulsion mediated by 
exchange. As a result, such a domain will not decay un- 
til its size reaches a microscopic scale (presumably on the 
order of a lattice constant), at which point the phase in- 
crease from to 2ir can be repaired through a phase slip. 
These nontrivial domains are the one-dimensional analog 
of the skyrmion^ a texture with a nonzero 0(3) winding 
number. Skyrmions may play a role in the magnetization 
reversal of thin films with high anisotropy^ 



APPENDIX A: PERIODIC STRUCTURE OF 
WIDE STRIPES 

When walls between uniform regions are well separated 
we can treat them as nearly free defects interacting with 
one another mostly via their stray magnetic fields. The 
stray fields of other walls are thus substantially weak and 
can be neglected when determining the structure of a 
well-isolated wall. Furthermore, its own stray field can 
be neglected as well: the energy associated with the stray 
field is 0(t 2 ), whereas all the other energies are 0(t). 
Finally, when the anisotropy highly favors out-of-plane 
magnetization, the characteristic width of the domain 
walls is small compared to that of the stripes, so that the 
energy to be gained by deforming the domain walls in the 
presence of an external magnetic field is far outstripped 
by the energy to be gained by moving them. Hence the 
external field may be neglected as well when determining 
the domain wall structure. 

It is convenient to use dimensionless variables for the 
effective anisotropy, field, coordinate, and wavenumber: 



K ~ h 
k = — , h = — , 



KQ 



k = 



kX 



In terms of these variables, the energy 
expressed as 



E 



local 



f 



fi Q M 2 L y t 2 4 



= - dx 



1 fd0_ 

2 \di 



cos 

2 



and 



E, 



stray 



fi Q M 2 Lyt 2 



' r ^\~k\\m z (k)\ 2 , 



(Al) 
2) can be 

hcos9 
(A2) 

(A3) 



where we use the Fourier transform of the out-of-plane 
magnetization 



(fc) = / d£e lte cos( 



(A4) 



We neglect for the moment the effects of the external and 
stray magnetic field in order to find the internal struc- 
ture of a domain wall. A domain wall that interpolates 
between cos 9 = — 1 and cos 9 = 1 and minimizes the sum 
of exchange and effective anisotropy terms in the energy 
above obeys: 



9' 2 = Ksin 2 i 



(A5) 
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The derivative here is with respect to x. This is solved by 
m z = cos 9 = ± tanh(v / ^x). The internal energy of each 
such wall is /j,oM 2 t 2 L y ^/~K/2, not counting the interaction 
energy due to the stray field. 

When calculating the stray field energy, it is easier to 
work with the i-derivative of the magnetization, and use 
the form: 



E, 



stray 



HoM 2 L y t 2 



dk |m' z (fe)| 2 



2tt 



Ifcl 



(A6) 
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If the walls have small spatial extent relative to the dis- 
tance between them, a periodic structure can be approx- 
imated as a sum of alternating upward and downward 
kinks. As such, we describe a periodic structure with 
period I — IX/ \/kq and upward length w — wX/^/kq in 
each period by the variational solution 

f d 

m' z (x) = / du A(w J , x — u) — (tanh(Vku)) , (A7) 

where 

00 

A(w,l,u)= ^ 8{u-nl)- 5{u-nl-w). (A8) 

n= — 00 



and 



ly/k sm(Trw /I) = it 2 exp 



(A15) 



Note that Eq. (|A15[) justifies the approximation that 
l-\Zhsxn.(7rw/l) 3> 1, since is large in this region of 
the phase diagram. 

These equations can be solved to give: 



w 



2 arccos(— h/ho) 



rr\/h$ - h 2 



and 



(A16) 



By substituting (|A7[) into (|A6|) . we arrive at 



E. 



stray 



167T 3 ^ . sin 2 (njw/ T) 
~W ^ 3 sinh 2 (tt 2 JfM) 



4(2^//) 3 d 
k dp 



-Ipjn 



j=0 n=l 



(A9) 



where p = tt 2 /\fkl and V — tL x L y is the volume of the 
film. 

Summing over j leads to an expression that can be 
approximated by an integral over n if p <C 1. We thus 
obtain 



E. 



stray 



PoM 2 Vk 



4 
Id 



ln(l + /) + 



2/ 
1 + / 



where 



f(wj) = 



l 2 k 



■ sm 



TTW 

T 



(A10) 



(All) 



The total internal energy of the kinks in this structure 
(including exchange and anisotropy) is 



PoM 2 Vk 



(A12) 



since there are two kinks in each period. The final contri- 
bution to the energy is the interaction with the applied 
magnetic field, described by the energy density 



E cx t. field 

PoAPVkq 



2hw 
I 



(A13) 



Minimizing the sum of these three terms with respect 
to / and w leads to the following expressions for the 
equilibrium period and upward width of the stripes [for 
l^/kam^w/l) 3> 1]: 



cos 



hir 2 



4a/I 



= ex P 



(A14) 



I = 



h 2 -h? 



(A17) 



with ho = — (A/7T 2 )Vkexp(— 7iVk/2). 

For h = the above equations reduce to w — 1/2 (i.e., 
no net magnetization), and a stable period of 



XI 



7T 2 A / ,T 

— 7= ex P 



(A18) 



Since the average magnetization (m z ) = 2w/l — 1, 
Eq. (|A14|) can be rewritten to give the magnetization 
curve for large k: 



2 . hir 2 
\m z ) — — arcsm — = exp 
tt \ 4Vk 




(A19) 



Further, note that as h—^—ho, the period I 
tends to infinity while the upward width stays finite 
[w — * (tt/4vk) exp(7rv^/2)]. For larger fields there is no 
stable periodic structure of this type. 

If, instead of a periodic structure, we are interested in 
the energy of a single soliton, we take I — ► L x 3> w in 
Eqs. (|A10p - (|A13p . so that f(w,l) » kw 2 /it 2 leading to 



4E 



t 2 L y p Q M 2 



iVk - 2hw 



const — 2hw 



In ^1 
In w 



2kw 2 



(A20) 



for an upward soliton of width w when kw 2 3> 1. 



APPENDIX B: SOLITONS AND THE 
EXCHANGE FORCE BETWEEN DOMAIN 
WALLS 



The energy of a single soliton given above by Eq. (|X20j) 
does not take into account the exchange interaction be- 
tween the domain walls bounding the soliton. This in- 
teraction will become important as the walls move closer 
together. 
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If we wish to find the effective force due to the exchange 
interaction, we must first find stable soliton solutions to 
the Lagrange equation associated with the energy (|A2|) . 
We will ignore the effects of the stray magnetic field (in 
the limit of large k with t <^ X/ y/K — \h\) . Since the 
force on the walls due to the magnetic field is known, the 
stable width of a soliton can be used to find the effective 
exchange force between the walls that must be acting to 
oppose the field. 

We will apply the boundary conditions 9 (±00) = 7r 
and 9' (±00) = 0, thus describing an upward soliton in a 
downwardly polarized background. We obtain from (|A2D 
that 



9' 2 1 

— = — iisin 9 — h(l 
2 2 v 



(Bl) 



where R = k/kq, h = h/no, and x = x^/ko/X, as in 
Eqs. (TMll-dAll). 

We expect different solutions to this equation for h > 
and h < 0. If h > 0, then the soliton is favored by the 
field, and the background is unfavored, held in place only 
by the out-of-plane anisotropy. In order to balance the 
force of the field, the exchange interaction in this case 
will attract the walls of the soliton to one another. 

The corresponding solution to Eq. (|B1[) at h > is 



cos 9 



2(1 - h/R) 



1 + (h/R) sinh 2 (fci) 



(B2) 



where k = yR—h, 

Note that cos 0(0) = 1 — 2h/R is the maximum value 
of cos 9 in this solution. Nowhere is there full upward 
polarization, so this solution describes a non-topological 
soliton, in which there is no net rotation of the magneti- 
zation between the ends of the system. 

If, on the other hand, h < 0, then the background is 
favored by the field and the soliton is not. The exchange 
force between the walls of the soliton must be repulsive in 
order to balance the force of the field squeezing the walls 
together. This is accomplished by a topological soliton 
solution, in which the magnetization rotates by a full 2tt 
between the ends of the system. The solution to Eq. (|B1[) 
in this case (h < 0) is 



2(1 - h/R) 



1 — (h/R) cosh (kx) 



(B3) 



Again, k = \J R — h. 

Note that in this case, cos 9(0) — 1 for any value of h. 
The 2-7T rotation of the magnetization forces the magneti- 
zation to point upward in the center of the soliton in any 
field. This is the reason for the repulsive exchange force. 
If the soliton is squeezed, the magnetization is forced to 
move from down to up and back in a shorter distance. 

If we set cos# = in Eq. ()B2|) to find the locations 
±Xo of the kinks bounding the non-topological soliton, 



we obtain 

h = 



1 + cosh 2 (xo'\/ R—h) 



4£exp {~2^x^j (B4) 



for large R. Since the force of the field is 2h acting to 
separate the kinks, the (rescaled) exchange force in the 
non-topological case is 



i^ cx = — 8«;exp ( — VR' 



(B5) 



Note that w = 2xq, since the kinks are at x = ±xq. 
Similarly, if we set cos 9 = in Eq. (|B3|) , we obtain 



- h = 



— 1 + sinh 2 (Jq yj R — h) 



4£exp (—2VRxq 



(B6) 

for large R and moderate width. Since the force of the 
field is now acting to compress the kinks, the (rescaled) 
exchange force in the case in which the soliton is topo- 
logical is F cx = -t-8£exp (^— v^Rw^j , a repulsive force. 

If the field is very high, the kinks in the topological 
case will be forced close together, so that the large Vkw 
approximation is no longer valid. The soliton will not 
collapse, however, until the width [obtained by solving 
Eq. JB6} for w = 2x ] 



wX 



2A 



arcsinh 



Vk+ \h\ 

is on the order of the lattice spacing 



1 



R 1 (B7) 



APPENDIX C: ENERGY OF STRIPES NEAR AN 
EDGE 

The energy of a two-dimensional thin-film system with 
magnetic field oriented normal to the plane of the sample 
and no bulk charge is 



E 



fioMH 



d z x 



A 2 |Vm 2 



2(1 
t [ d 2 k 



-ml) 2 
|fc||m z (k)| 2 . 



hm, 



(CI) 



Here m 2 (k) represents the two-dimensional Fourier 
transform of m z . In order to include the effects of a 
system edge, we will include a step function Q(y) as a 
factor in all the y-integrals that appear. For simplicity 
we discuss the case of zero applied field, h = 0, only. 

We use a trial solution with a single wave number and 
propagation direction. We will be comparing the energy 
of solutions with different values of the angle between the 
propagation direction and the system edge. We label this 
angle a. Our trial solution is 



asm (q • x — j3), 



(C2) 
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where a is the modulation amplitude, 
q= ((/cos a, gsina), and a ranges from for stripes 
running perpendicular to the edge to tt/2 for stripes 
running parallel to the edge. In order to evaluate the 
energy, we use the form 



e(y) 



dp ie ipy 
2tt p + ie 



(C3) 



for the step function in order to regularize the infrared 
divergence. Here L y = 1/e is the extent of the system 
in the y-direction. Where necessary, we also use the in- 
verse of the sample thickness as an ultraviolet cutoff in 
momentum space integrals. 



This integral requires both infrared and ultraviolet cut- 
offs. The inverse thickness 1/t serves to cut off the ul- 
traviolet divergence in the integral and e = 1/L y the in- 
frared. However, the difference E stlay (a, q) — E stray (0, q) 
is not sensitive to these cutoffs and can be evaluated in 
the limits t -> and 1/e = L y — > oo. Remarkably, the 
difference is also independent of the wavenumber q. By 
subtracting E stla y(0,q)/E o from Eq. (|C4|) and perform- 
ing the p integration, we obtain 



Astray q) ~ Astray (0, q) O 2 

E Q 8tt V ' 



(C6) 



where the function 



1. a /7r/2 

As long as stripes are not parallel to the edge, trial 
solutions with different values of (3 are related to one 
another by a translation in the x direction. Since the 
integration in Eq. (|C1[) is carried out over all x, the en- 
ergy is independent of the phase (3. Further, the total 
contribution of the local terms in the energy is the same 
for any orientation other than a = tt/2. The dependence 
of energy on the orientation then reflects the effect of the 
dipolar stray field alone. 

By inserting the trial solution (|C2j) into the stray field 
term in Eq. (|C1|) using the form (|C3j) for the step func- 
tion, we obtain directly that the contribution of the stray 
field to the energy is: 



Astray (a,g) _ _ f dp / y/g 2 + 2qp sj 

E 16 J 2tt V p* + e 



sin a + p 

2 



q 2 — 2qp sin a + p 2 

p2 _|_ £ 2 

where we introduce a characteristic energy scale 
E = h MH 2 L x . 



(C4) 



(C5) 



f(a) = sin a In (1 + sin a) + (1 — sin a) In (cos a) (C7) 

has a minimum at a = 0, when the stripes are normal to 
the edge. 

2. a = tt/2 

When stripes are parallel to the edge, the energy also 
depends on the phase (3. To leading order in t and e, 



E{n/2,g)-E(0,g) = 

E 8tt 



tt(k + X 2 q 2 ) sin(2/3) 
qt 



+ ln(2) - ln{qt) cos(2/3) .(C8) 

In addition to the energy of the stray field, this expression 
reflects the energy of anisotropy and exchange, both of 
which are sensitive to the positions of nodes relative to 
the edge. 

The overall energy dependence, including a dip at a = 
7r/2, is shown in Fig. [4] Depending on the parameters of 
the film, the global minimum may be cither at a = or at 
a = tt/2. However, even if the global minimum is at tt/2, 
the system may not easily find that configuration and 
remain in the metastable state with the stripes normal 
to the edge (a = 0). 
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